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Abstract 

The decay of the Standard Model Higgs boson in the intermediate-mass 
range into gluons is considered where special emphasis is put on the in- 
fluence of the leading electroweak corrections proportional to Gprn^. An 
effective Lagrangian approach is used where the top quark is integrated out. 
The evaluation of the coefficient function is performed using two different 
methods. The first one is concerned with the direct evaluation of the ver- 
tex diagrams and the second method is based on a low-energy theorem. 
In a first step the tools needed for the computation are provided namely 
the renormalization constants of the QCD Lagrangian are computed up 
to 0(a 2 s GFm%). Also the decoupling constants for the strong coupling con- 
stant a s and the light quark masses relating the quantities of the full theory 
to the corresponding quantities of the effective one are evaluated up to order 
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I. INTRODUCTION 



The Higgs boson is the only not yet discovered particle of the Standard Model of ele- 
mentary particle physics. Up to now only lower bounds on its mass of M# ;> 89.8 GeV [l 
could be derived from the lack of observation at LEP 1 and LEP 2. Through the virtual 
presence of the Higgs boson in loop diagrams it is also possible to set indirect limits on 
Mh with the help of the precision data collected at LEP, SLC and Tevatron. Currently 
they read M H = 84±l\ GeV with an upper limit of 280 GeV at 95 % C.L. [§. These num- 
bers suggest that a Higgs boson in the so-called intermediate-mass range, i.e. Mh ^ 2Mw 
where My/ is the mass of the W boson, is an attractive candidate. In this paper we will 
therefore consider such a Higgs boson and compute corrections to its gluonic decay rate. 

Comprehensive reviews concerning the properties of the Higgs boson are given in H§|. 
The dominant decay mode of an intermediate-mass Higgs boson is the one into bottom 
quarks. The QCD corrections are known up to 0(a^) |0-[7|. Concerning the electroweak 
theory, the full one-loop corrections are available ||. For the mixed electroweak/ QCD 
corrections only the leading terms of order a s Gpm\ || and a^Gpinf HlO] , |TTj| are evaluated 
at the two- and three-loop level, respectively. An other important decay mode is the one 
into gluons. However, this process is suppressed as compared to the fermionic channel 
as at lowest order it is mediated via a quark loop. It turns out that the leading order 



QCD corrections are quite large and amount to roughly 66 % [|T2|,|T3|]. Recently also 
the next-to-leading terms were evaluated which give a further enhancement of roughly 
20 % and thus increase the confidence to use perturbative QCD In J0| the leading 
electroweak corrections of order Gprnf were evaluated. In turned out that cancellations 
between different contributions take place and the final result is quite small. Nevertheless 
it is interesting to add an extra gluon and to evaluate the three-loop corrections of order 
a s Gi?mj, also in order to observe the behaviour of the perturbative series. Thus the aim 
of this paper is to consider next to QCD corrections also those terms arising from the 
electroweak theory which are enhanced by the top quark mass and proportional to Gpm\. 

For a Higgs boson in the intermediate mass range it makes sense to consider the limit 
M H <C ml and to construct in a first step an effective Lagrangian where the top quark 
is integrated out. Then the main task for the computation of the leading electroweak 
corrections is the evaluation of the effective coupling of the Higgs boson to gluons usually 
called C\. The operators appearing in the effective Lagrangian are only defined in the 
effective theory and therefore receive no corrections involving the top quark. Note that 
C\ enters not only into the decay rate but constitutes also a building block for the gluon 
fusion process which will be the dominant mechanism for the production of the Standard 
Model Higgs boson at the CERN Large Hadron Collider. The effective coupling to quarks 
is usually denoted by C 2 and was computed up to three loops in . 



It is well-known that the Appelquist-Carazzone decoupling theorem |16[ does not hold 
true in its naive sense if a renormalization scheme based on minimal subtraction is used. 
This means that the contribution from a heavy quark h with mass nih to a Green function 
of gluons and light quarks does in general not show the expected l/m h suppression. The 
standard solution to this problem is to do the decoupling "by hand" and to construct 
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an effective theory where the heavy quark is integrated out. The quantities relating 
the parameters, respectively, fields of the full theory to the corresponding quantities in 
the effective one are called decoupling constants. In [|I7| it was shown that there is a 
tight connection between the (renormalized) coefficient functions, C± and C2, and ( g , 
respectively, ( m which perform the decoupling for the strong coupling constant ot s and 
the light quark masses m q . Concerning pure QCD, ( g and ( m are known up to the three- 
loop level |T8|- pO| , [T7| . In this paper the leading electroweak corrections are considered and 
terms up to 0(a^Gpm^) are evaluated. The two-loop terms of order asGprn^ can be 
found in [[!]]]. 

The organization of the paper is as follows: In the next Section the notation is fixed 
and the theoretical framework developed in previous papers is reviewed. Section |T| is 
concerned with the computation of the renormalization constants for the QCD parameters 
and fields up to order a^Gpfnf. Although only the ones for the coupling constant a s and 
the light quark masses are needed we provide all renormalization constants of the QCD 
Lagrangian up to this order. In Section [TV] ( g and ( m are computed up to order o? s G F m\. 
Afterwards, in Section M, they are used in order to compute the coefficient functions 
C\ and C2. The result for C\ is compared with the direct evaluation of the triangle 
diagrams. C\ is then combined with the expectation values of the correlators formed by 
the corresponding operator in order to get a prediction for the gluonic decay rate of the 
Higgs boson. Finally, the conclusions are presented in Section [V] . 



II. THEORETICAL FRAMEWORK 

Let us in this section fix the notation and present the theoretical framework used for 
the calculation. The leading electroweak corrections are conveniently expressed in terms 
of the variable 

G F m 2 t 

where m t is the MS definition of the top quark mass which will be used throughout the 
paper. Corrections proportional to x t arise if in addition to the pure QCD Lagrangian 
also the couplings of the Higgs boson (h) and the neutral (x) and charged (0 ± ) Goldstone 
boson to the top quark are considered. 

For the evaluation of the decoupling constants it is necessary to know also the corre- 
sponding renormalization constants for the coupling and the masses 

9° s = H e Z g g s , m° q = Z mq m q , (2) 

up to the considered order. Here and in the following Zx denotes the renormalization 
constant in the MS scheme. Z g and Z mq will be computed in Section [ITT] together with 
the renormalization constants defined through 



^ - 1 = Z,{i - 1) , tf>° = J Ztlpq , 



G^ a = JZ Z G* c°' a = ^Z 3 c a , (3) 
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up to 0(a 2 x t ). g s = \J^o. s is the QCD gauge coupling, fi is the renormalization scale 
and D = 4 — 2e is the dimensionality of space time. m q is the MS mass of the light 
quark masses. G a is the gluon field, and c a is the Faddeev-Popov-ghost field. Colour 
indices for quark fields, ipg^ R , are suppressed for simplicity. However, it is necessary to 
distinguish between the renormalization mode of the right and the left handed quark field, 
ip L ' R = (1 ± 7 5 )^/2 as they are treated differently in the electroweak theory. For QCD 
we have, of course, Z% = Z R . The gauge parameter, £, is defined through the gluon 
propagator in lowest order, 

+ (4) 



q l + ie \ q 2 / 

The index "0" marks the bare quantities. Starting from the three-loop order, 0(a 2 Xt), 
the renormalization constant for the fermion wave function and mass, Z% and Z mq , 
depend on the quark species which is indicated by the additional index q. It represents 
one of the flavours u, d, s, c or b. Eqs. (|2]) and (|3]) also hold for q = t. If we refer to the 
four lightest quarks only the index q will be replaced by I. 

In analogy it is possible to write down the relations between the quantities in the 
effective (marked by a prime) and the full theory. Thereby we restrict ourselves to the 
case of the coupling constant and quark masses: 

9° s ' = cy s , < = «■ (5) 

Here, of course, q represents only one of the quarks u, d, s, c or b. It is more convenient 
to consider in a first step bare quantities and to perform the renormalization afterwards 
arriving at: 



<{y) = [^tQJ cl„{h) = ( 2 a s (n), (6) 

™>W = -^rCn q m q{^) = Cm q m q (fx). (7) 



7)1 



In the order considered in this paper the quantities in the effective theory are independent 
of the quark species. Thus, for the primed quantities the additional "q" is absent. A 
detailed derivation of the formulae for the computation of C° and was presented in fT?]. 



In turned out that only those diagrams where at least one top quark line is present have 
to be considered. They have to be evaluated for vanishing external momentum. The quite 
compact formulae for (® and Cm, read: 

Oh 



-o 



i + r G l(o,o) 



59 (l + iW))) Ji + n#(o) 



c o = hM ( 9 ) 

W 1 + S^(0) ' 



3 



where £y(p 2 ) and Es(p 2 ) are the vector and scalar components of the light-quark self- 
energy, defined through S(p) = ^(Ey(p 2 ) + 7 5 E^(p 2 )) + rriqUsij) 2 )- Note, however, that 
the axial- vector part Ha{p 2 ) does not enter explicitly into our analysis of ( mq as only 
leading order corrections in xt are considered [11]. The dependence of E(p) on "q" is sup- 
pressed. Il G (p 2 ) and n c (p 2 ) are the gluon and ghost vacuum polarizations and T G5c (p, q) 
is obtained from the one-particle irreducible diagrams contributing to the Gcc green func- 



tion p7| . The superscript "h" indicates that only the hard part of the respective quantities 
needs to be computed, i.e. only the diagrams involving the heavy quark contribute. The 
remaining paper is concerned with corrections of 0(x t ). Therefore in the following the 
full theory still contains the top quark, i.e. rif = 6 is the number of active flavours, and 
in the effective one the top quark is integrated out (rif = 5). Note that in contrast to 
the renormalization constants of Eqs. (0) and (^) the decoupling constants also receive 
contributions from the finite part of the loop integrals. 

In this paper the hadronic decay of a scalar Higgs boson in the so-called intermediate- 
mass range, i.e. M H <^ M w is considered. This process is affected by the virtual presence 
of the heavy top quark. Therefore it is promising to construct an effective Lagrangian 
which describes the coupling of the Higgs boson to (light) quarks and gluons. This has 
already been done in some detail in preceding works |2^ , p!2|j23| , pT|| . Hence only a brief 
sketch of the main steps and a collection of the relevant formulae is given. The starting 
point is the Yukawa Lagrange density describing the coupling of the H boson to quarks, 

£ y = -^ E ( 10 ) 

q£{u,d,s,c,b,t} 

where the sum runs over all quark flavours. In the limit m t — > oo Eq. (|10[) can be written 
as a sum over five operators p2[ , f23f formed by light degrees of freedom accompanied by 
coefficient functions containing the residual dependence on the top quark: 

v i=i 

It turns out that only two of the operators, in the following called 0[ and 0' 2 , contribute 
to physical processes. Expressed in terms of bare fields they read: 

0[ = (G^) 2 , O' = E «<■ (12) 

q£{u,d,s,c,b} 

The renormalized versions of 0[ and 0' 2 and the corresponding coefficient functions are 
given by: 




l + 2(a^/^)lnZ 



a 



(13) 



A(a' s d/da' s )\nZ' m ^ ^ 
C2 * " l + 2{ai s d/da' s )\nZ> Cl + °^ ■ (M) 
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Of course, the leading m t corrections considered in this paper only influence the coefficient 
functions as the operators are defined in the effective theory where no top quark is present. 

The factor H°/v° receives a finite universal renormalization and can be written in the 
form 



2^G X ' 2 H (l + S u 



(15) 



In [|K]] 5 U was evaluated up to 0(a 2 Xt). In our analysis only the 0(a s xt) terms enter. 
They are given by @: 



5 U = %t 



,(6) 



7 a 
2 + ~r; 



(16) 



( is Riemann's zeta function, with the value ((2) = ir 2 /6. 

Once the effective Lagrangian is at hand the evaluation of the hadronic decay rate splits 
into two parts, namely the computation of the imaginary part of the vacuum expectation 
values formed by the operators and the calculation of the coefficient functions. The 
correlators can be taken over form earlier works p4|,[7| as only pure QCD corrections are 
involved. In [0,0 two methods have been used for the computation of C%. The first one 

in the limit of a heavy 



relies on the direct evaluation of the triangle diagrams (see Fig. 
top quark. An expansion in the external momenta has to be performed up to linear order 
and the transversal structure has to be projected out. The second method is based on a 
low-energy theorem (LET) relating the coefficient functions to the decoupling constants. 



In []17j the following compact formulae were derived 



C x 



1 m t d 

2 dm 2 



InC 



C-2 



1 + 2 



mid 



dm 2 



InCn 



(17) 



which allow for a powerful check. Concerning the x t corrections it should be mentioned 
that the derivatives do not act on the overall factor m 2 . In this paper we use both methods 
in order to compute G\. Concerning 0%^ we use Eq. ( |TTD in order to reproduce the results 
given in [p] . 



III. RENORMALIZATION GROUP FUNCTIONS UP TO 0{a%X T ) 

In this section the renormalization constants for different parameters and fields of 
the QCD Lagrangian are computed in the MS scheme using dimensional regularization. 
Besides QCD corrections also the leading electroweak terms proportional to m\ are taken 
into account. Thus, in principle an index "(6)" indicating the number of active flavours 
would be necessary which is, however, omitted in this section. As only pole parts in D — 4 
have to be evaluated it is possible to reduce the calculation to massless propagator type 
diagrams, where the scale is given by the external momentum. In the case of pure QCD 
such a strategy is common practice, but also for the corrections of 0(x t ) this procedure 
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a 





LPJ)J)J)J)J)Jb 



FIG. 1. Feynman diagrams contribution to the coefficient function C\. The internal dashed 
line either represents the Higgs boson (h) or the neutral (x) or charged (^ ± ) Goldstone boson. 
The external dashed line corresponds to the Higgs boson. 



is applicable. Even in the presence of a heavy top quark the needful factors of m t either 
arise form the coupling of the scalar particles to the top quark or from the expansion of 
the internal top quark propagators at most up to linear order. 

The definition of the renormalization constants is given in Eqs. (0) and (^). In the 
following the computation of Z 3} Z 3 , Z\, and Z mq is presented. The renormalization 
constants both for the quark-gluon, the three- and four-gluon vertex can be obtained by 
the use of Slavnov- Taylor identities which relate them to the above five constants : 

* = = * = US) 

Zj 3 Zj 3 Zj 3 

For convenience we list in the following also the pure QCD results up to order a 2 . They 
were first obtained in | ]25|| . Note that from the results for the renormalization constants 
the corresponding anomalous dimensions can be computed. 

From Eqs. (|3]) it can be seen that for the computation of Z 3 corrections to the gluon 
propagator have to be considered. Some sample diagrams contributing at 0(a 2 x t ) are 
pictured in Fig. where the dashed lines represents either the h, x or ^ boson. The 
electroweak corrections arise for the first time at two loops. At three-loop order already 
224 diagrams contribute. Z% is obtained from the recursive solution of the equation 

Z 3 = l-K £ (u G (q 2 )Z 3 ) , (19) 

where Uq is the transversal part of the gluon polarization function defined through 

n^(g) = (-^V + g^)n G (g 2 ). (20) 



The operator K £ extracts the pole parts in e. 

After projecting out Il G (g 2 ) one ends up with purely massless diagrams. In this case no 
expansion in mt is required as the factor m 2 is provided by the coupling of the scalar bosons 
to the top quark. For the generation of the diagrams QGRAF is used. The output 
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FIG. 2. Feynman diagrams contribution to Z% and The dashed line either represents 
the Higgs boson (/i) or the neutral (x) or charged {(j)^) Goldstone boson. 



is then transformed to the package MINCER [27| which is written in FORM [25] and 
can deal with one-, two- and three-loop massless diagrams with one external momentum 
different from zero, so-called propagator-type diagrams. After the renormalization of the 
parameters a s ,m t and £ is performed for which only renormalization constants of lower 
order are needed Eq. (]TQ|) has to be solved recursively and one gets 



1 + 



a 



(6) i 



tt e 



/ 5 1 \ 1 1 /a (6) \ 



1 

+ - 



r 2 



25 



192 



V128 



384 



64 
1 



+ ^ + ttA )+ CaT (m + 1) - - — £ 



1 



+ x 4 T 



256" 128 
C A 



1 1 

4 + 8« 



48 24 

c f t (m + 1) ~ - c a t (m + i)^ 

£ V 4 A 48 



(21) 



where = c4 6) (/i). The superscript "(6)" indicates the number of active flavours. Ca = 
N c and Cp = {N 2 — 1)/(2N C ) are the Casimir operators of the adjoint and fundamental 
representation, respectively, and T = 1/2 is the trace normalization of the fundamental 
representation, m = rif — 1 = 5 is the number of light (massless) quark flavours. 

Special care has to be taken for the diagram pictured in Fig. 0(c). If the dashed line 
corresponds to the \ boson in each fermion line exactly one 75 matrix shows up and the 
naive treatment would fail. In turns out that the diagram is finite and thus gives no 
contribution to Z3. In the next section, however, an analogue diagram contributes to 
and a careful treatment is mandatory. 

In analogy to Eq. (|i~9| ) the renormalization constant for the ghost field is obtained 
form 



Z 3 = l-K e (u c (q 2 )Z 3 ) . (22) 
Corrections of order xt arise for the first time at three-loop level. The result reads: 
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1 

l + ^-C A 

n e 



1 J_ 

8 + 16^ 



a 



(6)' 



7T 



ci ( - 



1 

16 



256 



* + 512* 



+ c A T{m + i) 



32 



+ 



49 



1 



768 512 



a-C A T( ni + l) 



192 



+ x t C A T 



1 23 

8e2 + 967 



(23) 



The renormalization constant Z\ requires the computation of gluon-ghost vertex dia- 
grams, Tgcc{c[,p)- For simplicity we choose the gluon momentum to be zero and again end 
up with massless propagator-type integrals. This could in principle introduce unwanted 
infrared divergences. However, in it was shown that this is not the case and thus Z\ 
can be computed with the help of the formula 



(24) 



Also for Z\ the Xt corrections in principle appear for the first time at three-loop order. A 
closer look to the two-loop result shows, however, that those diagrams containing a closed 
fermion loop add up to zero. As the 0(a 2 x t ) terms come from the diagrams which contain 
a top quark loop accompanied with an additional exchange of a scalar particle we expect 
that Z\ gets no 0(a 2 Xt) corrections at all. This is verified by an explicit calculation. For 
convenience we list the pure QCD result as it is needed for the computation of Z g : 



Zt 



1 

l + ^-C A 

n e 

1 / 3 
+ - 

e 



4 + k 



a 



(6)' 



7T 



C 2 7 -t+±e 

A \e 2 \128 128^ 64^ 



7 



128 256 



1 

256 



(25) 



The charge renormalization constant, Z g , can be computed from the combination of 
Z 3 , Z 3 and Z\ with the result 



Zi 



Z 3 yZ^ 



1 + 



C 



it e 
121 



l 384 



11 11 

-C A — + T(ni + l)--x t T- 
24 v ; 6 2 



-C A T(n l + l)^- + T 2 (n l + V 2 ' 



a 



(6)' 



7T 



-C 2 A ^ + C A T( ni + l 



l 96 



+ x t T 



11 



48 v 24 

— + c f t (m + 1) — 

48 V ; 16 



-2 C A — + C F --T(ni + l)-)+-[ -G A - - C F - 



1 



(26) 



16 4 v ' 4/ e V 2 
As expected the £ dependence drops out which is an important check of the calculation. 
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From the fermion propagator the wave function renormalization constants Z^ R and 
the one for the mass, Z mq , can be computed. Due to the different coupling of the x boson 
to up- and down-type quarks, one gets an explicit dependence on the considered quark 
flavour which is indicated by an additional index. Although in Section |IV| only Z mq for 
light quarks is needed we consider for completeness also Z mt as the basic technique of the 
computation is very similar. The quark two-point function consists of three parts 



E(p) = i> E v (p 2 ) +7sSa(p 2 ) +m q E s (p 2 



(27) 



where Sy(p 2 ) and E^(p 2 ) correspond to the vector and axial-vector and £s(p 2 ) to the 
scalar part. Note that in the Standard Model the pseudo-scalar contribution to E(p) is 
zero. The renormalized fermion propagator can then be written as 



s-' 



(p) 



1-75 



-m. 



yL yR y 



1 - 



+ !±t*z? (1 



s s (p 2 



where S^(p 2 ) and £#(p 2 ) are given by 

S L (p 2 ) = Z v (p 2 )-Z A (p 2 ), E R (p 2 ) = 
From Eq. (p8| ) the following equations can be derived 

i-k £ (j: l (p 2 )z^ 

l-K £ {n R { v 2 )Z* 



E v (p 2 ) + Z A (p 2 



(2f 



(29) 



tR 



yL yR y 

Zjty Zjty Zj r , 



l + K E (E s (p 2 )^/z^zJ q Z mi 



(30) 



Their recursive solutions determine the wave function and mass renormalization constants. 

For the computation of Z^ R all masses appearing in the propagators may be set to 
zero and the factors mt are provided by the Yukawa couplings. In the case of the four light 
quarks only the diagrams where the gluons couple to the considered quark contribute and 
therefore the renormalization is left-right symmetric: 



1 + 



a 



(6) I 



it e 



-Of t£-t + 



a 



(er 



7T 



+ CaCf 

5 . 



1 

16 
1 

128^ 



64 s 64 s 



e)+c F T{ ni + r 



cl 



1 

e 
1 

16 



1 1 , 1 y 



64 

x t CpT 



32 

CaCf ( - 
1 1 



16" 32 
17 
64 



4s 



(31) 



Next to the pure QCD corrections the class of diagrams pictured in Fig. |3|(6) give rise to the 
0(a 2 Xt) corrections. In our approximation the bottom quark is effectively considered to be 
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FIG. 3. Feynman diagrams contribution to Z2 q , Z„ Lq and . The dashed line either 
represents the Higgs boson (h) or the neutral (%) or charged (cp^) Goldstone boson. 



massless which has the consequence that the bottom-specific corrections only contribute 



to ZJ: 
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+ C F T (m + l 
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16 
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CaCf 
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2e 
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24 
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(6)' 



7T 



64 



s 64 s 



+ -i) + C A C F 



64 

^ " tC(3) ) + CaC f 



/151 

V192 ~ 
31 

192 



64 



+ |C(3)) 



(32) 



with £(3) ~ 1.202 057. For the top quark altogether more diagrams have to be taken into 
account as — in contrast to the bottom case — already at one- and two-loop level the 
exchange of a Higgs and neutral Goldstone boson may occur. This is also the reason that 
both Z^ t and Z^ t get contributions from the top-specific diagrams. We get: 



7 L 
Z 2„ » 



(33) 



and 



Z£ + xA -- + 



a 



(6) 



c f 



7T 



2 " ^ + 



+ 



a 



(6)~ 



7T 
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Cf{ 16 + 2* 16* 



CaCf 



+ 



1 

+ - 

e 



32 s 64 s 



C F T (n, + 1) 



7_ 

24 



i-»+^(-| + |c(3)) + cr (ni + i)i 



(34) 



Note that the top-specific corrections of Eq. (|34D are exactly twice as large as the ones of 
Eq. ([§). 

As can be seen form Eq. (p0|) the scalar part can be used together with Z\ and Z^ q 
in order to compute Z mq . For it an expansion of £ (p) in m 3 up to linear order is necessary 
in order to be able to project out E5 and end up with massless two-point functions. Notice 
that here the factor m t may also origin from internal top quark propagators. The class 
of diagrams pictured in Fig. |(c) only contributes to E5. Special care has to be taken 
when the x boson is exchanged between the top quark loop and the light fermion line as 
then each fermion line contains exactly one 75. These diagrams, however, only develop 
an overall divergence giving rise to an simple \ je pole. Therefore we are allowed to adopt 
a prescription for 75 according to 't Hooft and Veltman [31] which was also used in 
and which is described in more detail in the next section. The result for q ^ b,t reads: 
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1 + — -C F 
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+ x t C F T 


■ 1 


2^2 




9 



11 



1 



- (C F — + C A C F — - C F T (m + 1) - 



32 



32 



+ 
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24 



3C(3)+2/3, -+3C(3) 



(35) 



The 0{a 2 s x t ) corrections arise from the diagrams shown in Fig. ^(6) and (c). 1 31 is the 
third component of the weak isospin, i.e. I 3 i = +1/2 for up-type quarks and I 3 i = —1/2 
for down-type quark flavours. The case of the bottom quark exhibits more structures and 
gives the result 
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(36) 



Finally for the case of the top quark we obtain 
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+C F T (m + 1) 
-C F T (m + 1) 



m 

32/ 

-) 

32/ 



+ 



(l|_9 c(3)) rtrv 



261 

128 ~ 
/239 

V128 



285 
128 



16 



C(3) 



(37) 



It is remarkable that except for the 1/s pole at 0(a 2 x t ) the coefficients of the structures 
x t (a^ / 7i) n / e m coincide in the expressions for Z mb and Z mt up to an overall sign. 



IV. DECOUPLING RELATIONS 

The computation of the decoupling constant relating a s in the full and the effective 
theory requires essentially three ingredients: the hard part of the gluon polarization func- 
tion, the one for the ghost polarization function and the one of the ghost-gluon vertex. 
At the one- loop order altogether only one diagram contributes to n<~ ; , namely the one 
containing a closed top quark loop which is obviously gauge invariant. This is also the 
case for the 0(a s x t ) result where inside the top loop an additional scalar particle is ex- 
changed. At two-loop level H® h and only receive pure QCD contributions. Actually, 
the individual diagrams contributing to give non-vanishing contributions, the sum, 
however, adds up to zero. At three-loop level there are 224, 14, respectively, 98 diagrams 
which have to be taken into account in order to compute of 11^, II^, respectively, r^ c 
up to order a 2 x t . Again, the separate diagrams contributing to r^g c add up to zero. 
gets non- vanishing corrections which have to be combined with For the generation 
of the diagrams the program QGRAF [EH] is used. The output is then transformed into 



a format suitable for the package MATAD []32| which is written in FORM [Ej| for the 
purpose to compute one-, two- and three- loop tadpole diagrams. 

A special treatment is necessary for the class of diagrams pictured in Fig. [|(c) where 
the dashed line represents the neutral CP odd Goldstone boson, x- Here, 7 5 occurs in 
two different fermion lines and the naive treatment would lead to a wrong result. Instead 
we follow the work of 't Hooft and Veltman [[HI] and write 7 5 in the form 



7 5 = ^W*7 M , (38) 

where >y^ upcT } is the anti-symmetrized product of four 7 matrices. The e tensor is pulled 
off from the analytical expression and an object containing eight indices is obtained: 

n£'£'(<7 2 ) = n Gil (g 2 )g^£'^<'4 ] + n G , 2 (g 2 )g [M g^^<'4 ] • (39) 
lie 1 and Ug 2 are functions of q 2 and may be extracted with the help of the projectors |3! 
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PL 



P%\xvpo (<?) 



24 (? 2 i^v4 ] - Hd^ui^ 



(q 2 ) 2 {D-l){D-2){D-3)(D-4) ' 
96 (-fgtftfgigt? + Dq^^gfg^ 



(1 



2\2 



(D — l)(D — 2)(D — 3)(D — 4) 



(40) 



Note that both Pi and P<i develop a 1/e pole which is a consequence of the fact that 
for D = 4 both structures appearing in Eq. (p9|) are linear dependent. This artificial 
pole, however, cancels in each diagram individually leading to a finite result. This was 
actually expected as there is no contribution to Z g from this class of diagrams at the order 
considered in this paper. 

Some comments concerning the renormalization are in order. The parameters in the 
lower order diagrams have to be replaced by the corresponding renormalized values. Con- 
cerning pure QCD, the counterterms for a s and mt have to be known up to order a s Xt- 
In principle also the QCD gauge parameter £ appears in the individual contributions Uq 1 , 
U® h and r^ c , however, only starting from three loops, which has no effect on the terms 
of order a 2 x t . The parameters which are present in the diagrams contributing to the 
0(a s x t ) results receive only contributions from pure QCD counterterms. 

In order to compute the renormalized quantity ( g also the renormalization constants 
Zf and Zf are needed. We choose to express the r.h.s. of Eq. (||) in terms of ay'. 
Therefore Z^ gets a dependence on m t through the substitution of whereas before 
only pure QCD terms were present. Z^ has to be known up to order a 2 Xt which was 
derived in the previous section. Finally we obtain for Q the following result: 



c 2 



2 , °i 6) M T 




TT 



+ riin 2 ^ 



mi 



25 




C F -^ + -lnAr 



+ x t 



C(2) + — C(3)-31n^ + -ln 



mi 



4 



,2 V_ 

9 

mf 



+ Ca ( + -C(2) - — C(3) + - In ^4 ) + T [5 + Z C (3) + i In r- 
1 4 8 SV ' 64 sv ; 4 m 2 t U 8 SV ; 9 m 2 t 



3 mf 



7 



C(3)T 



(41) 



where the contribution of the diagrams in Fig. 0(c) corresponds to the last entry in the 
last line of Eq. (f|l]). For convenience also the pure QCD result of 0(a 2 ) is listed. The 
corresponding three-loop terms can be found in |54]]I7H ■ 

The decoupling constants for the light masses, ( mq , requires the computation of the 
hard part of the fermion propagator. Implicitly this has already been done in [IT] where, 
however, the main focus was on the evaluation of C^, the effective coupling to light 
quarks. In this paper we will also list the result for ( m . As already mentioned above 
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Cm 9 depends on the considered quark flavour which can be seen by a look to the diagrams 
pictured in Fig. ^. Some words are in order in connections with the diagram in Fig. §(c) 
when the dashed line corresponds to the x boson. Actually this diagram is responsible 
for the difference between ( mu and ( md , as the xff coupling is proportional to the third 
component of the isospin. Furthermore, the treatment of 75 needs some care. As already 
mentioned in Section |T| the diagrams of this class have an overall divergence which is 
reflected in the 1/e pole contributing to Z mq whereas all subdiagrams are finite. Thus we 
are allowed to adopt the prescription for 75 described above. The decoupling constant for 
the u, d, s and c quark then reads (for simplicity we set N c = 3): 



c 



HI, 




36 



1 1 2 A* 
+ — In 2 



mi 



12 



mi 



>) + gc(3) 



T 2 2 

9C(4)--ln^ + 6C(3)ln^ 



6 



mi 



mi 



+2hi 

For the bottom quark one receives: 



51-^(3) + 9C (4)- in 



6C(3)ln£ I 



mi 



(42) 




5 3, u 2 a® 
- + - In At + — 



(M) 



mi 



7T 



16 A/ x 7 u 2 3 , 2 u z 
--4C(2) + -ln^ + -ln 2 £^ 



mi 



mi 



472933 



12096 

55 



6133 A/ . 905 A/ . 383 A/ , N 1251 19 
<(2) + — C(3) + — C(4) + — S 2 + —D 3 



168 
43 



C(2)--C(3) ln£,+ 



72 
K 2 



112 



72 



mi 



529 
48~ 



In 



2 V> 



+ 



29 
12 



In 



3 » 



C(2)-2C(3) + 



•- + -C(2) 
144 3 s V ; 



In 



-^ln 2 ^ 

2 m? 



2 



-^ln 3 ^ 



where we have used Cp 
4 



(43) 



4/3, Ca = 3 and T = 1/2. The constants 



S 2 
B 4 



9V3 



0.260 434 



6C(3) - H C(4 ) 



6 Cl 2 - 



7T 



-3.027009 



16Li 4 (- 



> 



4C(2)ln 2 2 + -ln 4 2 



-1.762 800 



(44) 



where ((4) = 7r 4 /90, CI2 is Clausen's function and Li 4 is the quadrilogarithm, occur in the 
evaluation of the three- loop master diagrams P5Vf?7|] . In |T7| the three-loop corrections of 
0{a z s ) were computed. The results of Eqs. (f|I|) and fl43D will be used in the next section 
in order to compute Ci x and C 2b - 
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V. HADRONIC HIGGS DECAY 



In this section we compute the coefficient function C\ using two different methods. 
The first one is concerned with the direct evaluation of the triangle diagrams (see Fig. |l|) 
connecting the Higgs boson to two gluons. In a second step C\ is evaluated with the help 
of the LET where the result of the decoupling constant ( g derived in the previous section 
is used. 

At three-loop level altogether 990 diagrams contribute to the order o? s xt- Some typical 
examples are pictured in Fig. [I] where the internal dashed line either represents the Higgs 
boson, h, the neutral Goldstone boson, Xi or the charged Goldstone boson, cf) ± . In the 
latter case also the bottom quark is present in the fermion loop. The external Higgs boson, 
however, only couples to top quarks as corrections proportional to x t are considered. The 
diagrams have to be expanded in both external momenta and the transversal structure 
which arises from the external gluons is projected out in order to end up with scalar 
integrals. 

Again, the packages QGRAF p6 and MATAD [3^] are used for the generation, re- 



spectively, the computation of the diagrams. A general gauge parameter, £, for QCD[] is 
used and the independence of the final result serves as a welcome check for the correctness 
of the result. There is again a class of diagrams which requires a special treatment due 
to the fact that 75 appears in two different fermion lines. A sample diagram is shown 
in Fig. |T](c) . Of course, this class is tightly connected to the one discussed in connection 
with ( g (see Fig. 0(c)) and we can adopt the handling for 75 developed in Section [TV]. 
From Eq. (f|l]) one can see that no \nm t terms are present which has the consequence 
that according to the low-energy theorem of Eq. ([H]) C\ gets no contribution from these 
diagrams. This is confirmed by the direct calculation of the vertex diagrams: there is no 
contribution from this class of graphs. 

After taking into account the counterterms needed to get the renormalized coefficient 
function (see Eq. (ITj) ) one arrives at: 



^ifWL.^L 3 c 5 T i j* 

7r 71 4 4 3 mi 



+ x t [C F 9 




2 ml 



(45) 



The 0{a s xt) terms can be found in |L5[] and the 0(a 2 s ) results were computed in fi~2|Ji! 
The QCD corrections of 0{a i s ) and 0{al) are also known [13,17], however, for simplicity 
they are not displayed in Eq. (|4"5j) . 

If we use the decoupling constant of Eq. ([|l|) and plug it into Eq. (|T7|) we obtain 
the identical result which serves as a non-trivial check. Note that the diagrams to be 
considered in both approaches are quite different. Furthermore, the method based on the 



In the considered limit the electroweak gauge parameters drop out trivially. 
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LET requires the renormalization constant Z g to be known at 0(a 2 x t ) whereas for the 
direct computation only the terms up to order a s x t are necessary. 

For completeness we also list the result for C 2 obtained from Eqs. (f42|), respectively, 
((43p and Eq. (|17|). For the light quarks we get 



C 2 



1 + 



7T 



11 



\ In + x t 

3 mi 



7 



12C(3) + 2/ 3i (2 + 12C(3)) 



(46) 



The coefficient function for the bottom quark reads: 



C 2b = C 2d + x t < 



7T 
,2 



-7 -6 In 



mj 



+ 



7T 



12169 110 w x 



43 89, (S 55, 2 fi 2 

2 m. 



^ 72 3 SV ; m? 



(47) 



In [|TTJ C 2; and C 2b are listed for general gauge group SU(N C ). 

Let us now have a look at the numerical consequences on the decay rate. Therefore 
also the imaginary part of the correlator ([O'JfC^]) is needed which can be found in \TA . 
Furthermore the universal corrections arising from 5 U (see Eq. ([15])) have to be taken into 
account. Inserting all building blocks into the equation 



M, 



1 + 5 U ) Cllmip'M]) 



and expanding up to the three-loop level leads to 
r (H - gg) 



r Born (H - gg) 



1 + x t + s v ; [23.750 - 1.667nj 



7T 



+ x t 38.837 + 2.000 In 



Ml 
m 2 t 



1.167nj 



+ 



af\M H ) 



7T 



370.196 - 47.186nj + 0.902n 2 + (2.375 + 0.667n,) In 



M 



11 



m, 



l + x t + 



(M; 



7T 



17.917 + x t 33.004 + 2.000 In 



+ 



2 r 



7T 



M 2 

156.808 + 5.708 In — ^ 

mi 



(4f 



(49) 



with r Born (# -> (/(/) = G F M^/mKy/2 x (af)(M H )^) 2 . The renormalization scale /x 2 
is set to Mjj. The resulting logarithms In Mjj/m 2 are numerically small which makes 
a resummation not necessary [14,33|. In a similar way to the leading electroweak cor- 
rections [15] also at 0(a s xt) large cancellations between the universal and non-universal 
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terms take place. Actually, the "1" in front of the 0(x t ) term is composed of (7—6) and the 
"33.004" results from (131.504 — 98.500) where in both cases the first number corresponds 
to the universal corrections. Choosing m t = tu^Mh) = 173 GeV and Mh = 100 GeV 
leads to 

J^M_ + + (17 , 17 + 3o,n,«) + 15 o,5o f») " 

« 1 + 0.0031 + 0.6559 + 0.0035 + 0.2017 . (50) 

If one compares the 0(xt) corrections to the pure QCD terms of the same loop order they 
are clearly negligible. The three-loop corrections of 0(a s xt) amount roughly 1 % of order 
a 2 s terms. It is, however, interesting to note that they have the same sign and that they 
are of the same order of magnitude as the two- loop corrections of order x t . 



VI. CONCLUSIONS 



In this paper the gluonic decay of a scalar Higgs boson in the intermediate mass range 
is considered. The pure QCD corrections up to 0{a A s ) were computed recently |14| and it 



turned out that they are quite sizeable. On the other hand the leading electroweak cor- 
rections are small. The main focus of the present paper was the evaluation of additional 
QCD corrections to the O(xt) term. Although the new corrections are small as compared 
to the QCD terms they are quite sizeable as compared to the leading electroweak correc- 
tions. It seems that once QCD is switched on and its full structure is available, i.e., the 
gluon self-interaction is at work, large corrections can be expected. This is also motivated 
by the comparison of the leading order QCD corrections to T(H — > gg) and T(H — > 77). 
In the latter case only the Abelian part of the QCD enters and the correction factor only 
amounts to 1 — 2 % . 

In this paper the top quark is considered to be much heavier than the other quarks 
and the Higgs boson. However, the generalization of the analysis to any heavy quark 
which fulfills these conditions is straightforward. 
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